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Morphometric variability in round goby Neogobius 
melanostomus (Perciformes: Gobiidae) from the Sea of Azov 
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A total of 38 morphometric characters of the round goby Neogobius melanostomus melanosto- 
mus (Pallas, 1814) were studied in its native range in the Sea of Azov. The aim was to assess 
joint effect of sexual dimorphism and size variability on overall variability within groups of 
samples (populations) using appropriate methods of traditional statistical analysis (one-di- 
mensional and multivariate statistics). Sex and size-dependent variability was studied based 
on model samples of males and females of different size. Most of the studied morphometric 
characters of round goby from the Sea of Azov demonstrated statistically significant sex- and 
size-dependent in-group variability. The pattern of the variability suggests that, for a com- 
parison of round goby from different sea regions and between populations, separate samples 
of males and females within the range 9—13 cm SL should be examined in order to minimise 
the effect of the size and sex factors. The approach of searching for the most informative size 
range could be useful not only for further studies of infraspecific variation but for comparisons 
between morphologically close gobiin species. 


С целью оценки совместного влияния полового диморфизма и размеров рыб, используя 
соответствующие методы статистического анализа (одномерные и многомерные) изуче- 
на морфологическая изменчивость бычка-кругляка Neogobius melanostomus melanostomus 
(Pallas, 1814) из Азовского моря — нативного ареала вида. Всего изучено 38 пластиче- 
ских признаков самок и самцов разных размеров. Подавляющее большинство изучен- 
ных морфологических признаков бычка-кругляка Азовского моря демонстрируют ста- 
тистически значимую половую и размерную изменчивость. Характер морфологической 
изменчивости показывает, что при сравнительных исследованиях, с целью минимизации 
влияния размеров и пола рыб, анализируемые выборки должны включать рыб одного 
пола в интервале длин (SL) 9—13 см. Аналогичный подход, основанный на выявлении 
наиболее информативного размерного интервала, может быть использован не только 
для изучения внутривидовой изменчивости, но и в таксономических исследованиях 
близких видов бычковых рыб. 
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INTRODUCTION 


The native range of the round goby Ne- 
ogobius melanostomus melanostomus (Pal- 
las, 1814) includes basins of the Black and 
Marmara seas and the Sea of Azov where 
it inhabited mostly lower reaches of rivers 
and the whole brackish Sea of Azov. In the 
latter, the round goby is one of the most 
abundant and commercially valuable spe- 
cies (up to a third of annual fish harvest in 
2008—2015). Neogobius melanostomus af- 
finis (Eichwald, 1831) distributed in the 
Caspian Sea basin, mostly inhabits marine 
and estuarine habitats. 

Round goby are spread across the sea 
from almost freshwater areas with a salinity 
of about 0.5% (eastern Gulf of Taganrog) 
to over 15%о (near the Kerch Strait). The 
most dense aggregations of round goby in 
summer and autumn, when they intensively 
forage, occur in areas with high biomass of 
zoobenthos with a predominance of bivalve 
molluscs in the western part of the sea. 
Round goby do not perform large-scale mi- 
grations and only move to coastal areas in 
spring and away from the coast to slightly 
deeper areas in autumn. 

In the second half of the last century, 
the round goby began to expand the area 
of distribution upstream the rivers Danube, 
Dniester, Southern Bug, Dnieper, Don, and 
Volga (Smirnov, 2001; Pinchuk et al., 2003; 
Bogutskaya et al., 2004). In 1990, with the 
ballast water of ships the round goby was 
introduced into the Baltic Sea (Skóra, Sto- 
larski, 1993), and began to actively spread, 
first, in the Baltic Sea (Ojaveer, 2006) and 
then in the coastal zone of the North Sea. 
The round goby was introduced with the 
ballast water of ships into the North Ameri- 
can continent where it was first collected in 
1990 and then became successfully estab- 
lished and currently occurs in the basin of 
all the Great Lakes (Jude et al., 1992; Cor- 
kum et al., 2004; Kornis et al., 2012). So 
far, the round goby has become an object of 
intensive research due to a significant ex- 
pansion of its nonindigenous range and the 
impact on native aquatic ecosystems. 
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Morphology of the round goby in its na- 
tive range is not sufficiently studied; in par- 
ticular, age-and-size variability and sexual 
dimorphism are poorly known which is crit- 
ical to ensure comparability across samples 
and populations. Some aspects of morpho- 
logical variability and sex-dependent mor- 
phometric differences in round goby have 
been studied earlier. Thus, it was shown that 
some morphometric characters correlated 
with body size and sex in the round goby 
from the Dnieper-Bug estuary (Bil'ko, 1971; 
Bil'ko, Vybornaya, 1972). Smirnov (1986) 
reported the presence of sexual dimorphism 
in this species in the Molochny Liman estu- 
ary and the Danubian coast of the Black Sea 
but found no difference between males and 
females in the Sasyk Lagoon. 

The absence of sex-dependent differ- 
ences was also reported for the round goby 
outside its native range, e.g., the Kuibyshev 
and Saratov reservoirs (Shemonaev, 2006) 
and the middle reaches of the Danube 
River (Simonovié et al., 2001). These con- 
clusions were based solely on the results 
of statistical analyses using Student's test 
but without providing two mandatory con- 
ditions for its applicability (Lakin, 1980; 
Orlov, 2002): 1) the use of Student's test 
assumes a normal distribution in samples; 
this assumption can only be violated with 
large samples, and 2) compared samples 
need to have the same variance. However, 
the publications mentioned above did not 
discuss the question of applicability of Stu- 
dent's t-test, even when very small samples 
were compared. 

Data on sexual dimorphism and size var- 
iability of morphometric characters in the 
round goby from the Sea of Azov were first 
published only recently. Zabroda (2009) 
and Zabroda & Diripasko (2009) reported 
on size-dependent variability in males and 
compared morphometric characters in 
males from two different localities in view 
of the revealed features of size variability. 
These data demonstrated further needs for 
a comprehensive and integrated approach 
to assess joint effect of sexual dimorphism 
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and size variability on overall variability 
within groups of samples (populations) and 
variation between groups using appropri- 
ate methods of statistical analysis. So, the 
aim of this study was to analyse sex and size 
variability based on model samples of round 
goby males and females of different size. 


MATERIAL AND METHODS 


The material, containing 50 females 
and 65 males, was collected in early spring 
(when matured ovaries do not affect female 
body shape) in 2009 from catches of gob- 
ies done by mechanised dredges in North- 
Western region (at a location at 46"10^N, 
35°20’E) of the Sea of Azov. 

As round goby males attain larger maxi- 
mum size than that in females, the statisti- 
cal analyses were performed on both males 
and females within one and the same stand- 
ard length (SL) range, 5-15 cm. Males of 
the entire SL range (5-18 cm) were also 
analysed. 

A total of 40 measurement were done 
(Fig. 1) using a calliper to a nearest 0.1 mm. 
For the subsequent statistical processing, 
relative morphometric characters as per- 
centage of SL and head length (HL) were 
used. 

For statistical processing of data, meth- 
ods of one-dimensional and multivariate 
statistics were used. 

In order to be able to compare the re- 
sults of different methods, both parametric 
and nonparametric statistical tests were ap- 
plied. As the first step, a Shapiro-Wilk nor- 
mality test was performed to check if the 
variables (morphometric characters) follow 
a normal distribution, a mandatory condi- 
tion to be met for using parametric tests. 

Nonparametric Mann-Whitney test 
(U-test) and parametric Student's test 
(£-test) were used to infer about statistical 
significance of differences between the mor- 
phometric characters. 

Statistical correlation between the char- 
acters and SL was measured using two coef- 
ficients: 1. Pearson's correlation coefficient 
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for normally distributed variables and 2. 
Spearman's rank correlation coefficient for 
unknown or non-normal distributions be- 
cause it is less sensitive to non-normality 
(Sachs, 1976; Rebrova, 2002). Correlation 
coefficient (Pearson's r) absolute values 
from 0.50 to 0.75 indicate moderate to good 
correlation, and r values from 0.75 to 1 indi- 
cate very good to excellent correlation be- 
tween variables (Dawson & Trapp, 2004). 

To evaluate statistical dependence of 
morphometric characters on SL and sex, the 
following methods were also used: 

1. Multivariate analysis of variance 
(MANOVA). Despite the fact that this 
method is, in fact, parametric and does 
not have any immediate nonparametric 
analogues, it is applicable as being not par- 
ticularly robust to departures from multi- 
variate normality (Anderson, 2001). The 
MANOVA analysis was preceded by testing 
for distribution normality of either of mor- 
phometric variables and these data were 
taken into account when interpreting the 
results. 

2. Univariate analysis of variance 
(ANOVA) which analyses the differences 
among group means and their associated 
procedures such as descriptive statis- 
tics and correlations between variables 
among and between groups (in our case, 
size groups with a 2 cm SL interval). Since 
this method of test statistics are applicable 
for parametric statistical hypotheses, the 
Kruskal-Wallis test, which is a nonpara- 
metric analogue of ANOVA, was also used. 

For the purpose of classification of sam- 
ples, cluster analyses (using the Complete 
Linkage method and the Euclidean dis- 
tance as a proximity/similarity measure) 
were implemented using Microsoft Excel 
and Statistica 6.0 software. 

Abbreviations: n, number of specimens; 
SD, standard deviation; r, Pearson's corre- 
lation coefficient: R, Spearman's rank cor- 
relation coefficient; U-test, Mann—Whit- 
ney’s nonparametric test; t-test, Student's 
test. For abbreviations of measurements see 
Fig. 1. 
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Fig. 1. Scheme of morphometric measurements round goby: a) side view, b) top view, c) a bottom 
view. 1—2, standard length (SL); 1—3, total length (TL); 1—4, predorsal distance (aD); 1—9, prepec- 
toral distance (aP); 1—10, prepelvic distance (aV); 1—11, preanal distance (aA); 1—26, head length 
(HL); 1—29, preorbital distance (ao); 1—31, upper jaw length (1m); 1—32, lower jaw length (Imd); 
2—3, caudal fin length (1C); 2—8, postdorsal distance (pD); 2-13, caudal peduncle length (pl); 4—5, 
maximum body depth (Н); 4—14, length of first dorsal fin base (11D); 6-7, minimum body depth 
(caudal peduncle depth) (h); 9—10, pectoral-pelvic distance 1 (from upper end of pectoral fin base) 
(Рет) У); 9—12, width of pectoral fin base (iP); 9—24, pectoral fin length (1P); 10—11, pelvic-anal 
distance (V—A); 10—12, pectoral-pelvic distance 2 (from lower end of pectoral fin base) (Pian V); 
10—25, pelvic fin length (ТУ); 11—21, length of anal fin base (1A); 15—16, first dorsal fin depth (h1D); 
17—18, length of second dorsal fin base (12D); 19—20, second dorsal fin depth (h2D); 22—23, anal fin 
depth (hA); 26-30, postorbital distance (ор); 27—28, head depth at паре (hcz); 29-30, horizontal 
eye diameter (о); 33—34, distance between eye and corner of mouth (or); 35—36, cheek depth (hop); 
37—38, head depth through middle of eye (hco); 39—40, maximum body width (iH); 41—42, mini- 
mum body width (caudal peduncle width) (ih); 43—44, head width (ic); 45—46, interorbital distance 
(io); 47—48, mouth width (ir); 49—50, width of pelvic fin base (iV); 51—52, isthmus width (ist). 
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RESULTS 


Data of statistical analyses applied for 
two samples, males (49 specimens) and fe- 
males (50 specimens), both groups includ- 
ing fish 5—15 cm long, are shown in Tables 
1 and 2. Results of the analysis of size-de- 
pendent morphometric variability of males 
within the entire range of SL (5—18 cm) are 
given in Table 3. 

The results of the Shapiro-Wilk test 
showed that most of the examined charac- 
ters are normally distributed; the charac- 
ters demonstrating some departures from 
the normal distribution either in males or 
in females are marked with # in Tables 1—3. 

The degrees of correlation estimated by 
the Pearson's correlation coefficient and 
Spearman's rank coefficient are very simi- 
lar. The exceptions were two characters in 
males (Table 3): minimum body width and 
second dorsal fin depth. 

For the minimum body width index 
which departs from the normal distribution 
in males, the Pearson's coefficient is statisti- 
cally significant (р=0.05) while the Spear- 
man's coefficient is not. For the second 
dorsal fin depth index which is normally 
distributed, the results are opposite. 

For most of the examined characters, a 
statistically significant correlation of mod- 
erate to good power (r absolute values 0.5— 
0.75) was revealed while some characters 
are strongly (r > 0.75) correlated with SL 
(Tables 1—3). The direction of correlation 
(either positive or negative) of all characters 
with statistically significant correlation co- 
efficients is the same in the males and in the 
females. However, characters strongly cor- 
related with SL are not entirely the same in 
the females and in the males; in the females, 
ten of them are strongly correlated with SL: 
five describing the body shape (?6 SL) and 
five describing the head shape (?6 HL), but 
in the males (similarin both groups, 5—15 cm 
and 5—18 cm), eight characters are strongly 
correlated with SL: one describing the body 
shape and seven describing the head shape 
(Tables 1, 3). The results of the variance 
analyses (ANOVA and Kruskal-Wallis) of 
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round goby males including the large-sized 
specimens were similar except for a single 
character (predorsal distance) (Table 3). 

The evaluation of differences in charac- 
ters depending on sex of fish using the two 
tests resulted in similar values (Table 1). 
Only two characters (maximum body depth 
and lower distance between pectoral and 
pelvic fins) were found to be different with 
regard to the level of significance obtained 
from the U-test and the t-test: the t-test re- 
vealed no statistically significant differenc- 
es (p > 0.05) while the U-test showed sta- 
tistically significant differences (p x 0.05) 
(Table 1). In total, the males and females 
differ by 21 from 38 examined characters. 

Thus, most of the studied morphometric 
characters of round goby in the examined 
samples from the north-western region of 
the Sea of Azov demonstrate statistically 
significant sex- and size-dependent in- 
group variability. 

The results of the MANOVA (Table 2) 
showed that the vast majority of the mor- 
phometric characters depend on both sex 
and SL. Only three of them (pD, 12D, and 
op) out of 38 analysed, are not confirmed to 
besignificantly dependent on the size of fish 
and/or sex (Table 2). 

The variability dependent on both size 
and sex of fish based on the entire set of ex- 
amined morphometric characters is shown 
by a dendrogramme resulting from a clus- 
ter analysis (Fig. 2). A separate branch of 
the graph shows a group of juvenile round 
goby of 5—7 cm SL of both sexes. A distinct 
cluster is formed by large males with SL 
exceeding 13 cm (three size groups: 13—15, 
15-17, and 17-19 cm). This cluster does 
not include females of 13-15 cm SL that 
means that they are more similar smaller in- 
dividuals of both sexes than to males of the 
same size. Also, the females and males in the 
size group 7—9 cm are much more similar to 
each other than to individuals of their re- 
spective sex from the adjacent size groups. 
The most similar are the males from the size 
groups 9—11 cm and 11—13 cm and the fe- 
males from the same size groups (Fig. 2). 
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Table 2. Morphometric characters of round goby from north-western Sea of Azov and results of cor- 
relation and MANOVA analyses to assess dependence on sex and SL factors. 























Males and females Correlation MANOVA analysis, 
(n = 99) analysis p actual 
Character 
Mean SD Min Max r R Sex SL vie 
and SL 

SL, cm 10.14 2.747 5.1 14.8 

% SL 

H 23.24 1.721 1804 27.97 | 0.621 0.588 0.001 0.000 0.419 
h 11.55 0.731 9.45 13.85 | 0.612 0.619 | 0.001 0.000 0.707 
iH 18.93 1489 15.19 21.90 | 0.598 0.574 | 0.000 0.000 0.351 
ih #5.14 0.441 4.36 6.33 | -0.326  -0.335 | 0.007 0.000 0.023 
ар 32.78 1158 2938 35.61 | 0.362 0337 | 0.237 0.000 0.144 
рр 19.61 1.143 16.85 22.22 | -0.173 -0.159 | 0.079 0.398 0.689 
aP 29.31 0.999 2682 31.73 | 0.409 0400 | 0.018 0.000 0.068 
ау 29.81 1.306 26.67 32.65 | 0.507 0481 | 0.000 0.001 0.104 
аА 55.74 1.902 51.18 59.66 | 0.384 0362 | 0.809 0.018 0.303 
Po)-V #1638 1493 1241 19.08 | 0.792 0.780 | 0.617 0.000 0.433 
PV 5.91 0.691 3.89 7.48 0.555 0.527 | 0.000 0.000 0.003 
V-A 27.32 1.939 22.04 32.03 | 0.295 0.260 | 0.068 0.005 0.084 
р1 21.02 1.195 18.80 24.80 | -0.246 -0.238 | 0.593 0.025 0.573 
Ир 17.13 1419 13.54 19.77 | 0454 0422 0.002 0.000 0.417 
һр 15.16 1.102 12.24 17.72 | -0.004 -0.017 | 0.000 0.520 0.000 
12р 32.80 1.356 30.00 36.83 | 0.256 0.274 | 0.516 0.111 0.499 
h2D 15.77 1017 13.22 18.70 | -0.310 -0.281 | 0.683 0.016 0.002 
lA 25.21 1.283 22.15 2881 | 0.173 0.125 | 0.164 0.297 0.043 
hA 311.36 0.836 9.62 14.07 | -0.397 -0.344 | 0.002 0.000 0.950 
IP 30.29 1.499 26.56 3433 | -0.471 -0.504 | 0.503 0.000 0.128 
iP 12.24 0.815 10.00 1441 | 0.707 0.695 | 0.142 0.000 0.000 
IV #2359 1.243 20.00 27.22 | -0.549 -0.541 | 0.037 0.000 0.830 
iV 7.01 0.512 5.64 8.20 0.709 0702 | 0.007 0.000 0426 
IC 21.47 1.267 18.72 25.16 | -0.599 -0.588 | 0.246 0.000 0.145 
HL 27.33 1.070 25.00 30.19 | 0.269 0.290 | 0.000 0.010 0.023 
% HL 

hez 81.62 5.530 6438 92.62 | 0.726 0.718 | 0.002 0.000 0.085 
ic #90.44 7.575 65.75 101.46 | 0.799 0.796 | 0.031 0.000 0.586 
ao #31.24 2.593 23.65 38.17 | 0.626 0.644 | 0.341 0.000 0.758 
о 22.54 2094 16.95 27.63 | -0.794  -0.792 | 0.005 0.000 0.386 
ор 51.19 1.700 46.01 5479 | 0.293 0270 | 0.097 0.111 0.323 





© 2017 Zoological Institute, Russian Academy of Sciences, Zoosystematica Козяса 26(2): 392—405 


400 О.А. DiRIPASKO € Т.А. ZABRODA. MORPHOMETRIC VARIABILITY IN NEOGOBIUS MELANOSTOMUS 


Table 2. Continued. 

















Males and females Correlation MANOVA analysis, 
(n = 99) analysis p actual 
Character 
Mean SD Min Мах r R Sx SL 5% 
ка а к апа SL 
10 17.28 2672 9.87 23.36 | 0.835 0841 | 0.016 0.000 0.852 
Im #31.68 3.280 2424 38.17 | 0.762 0.767 0.001 0.000 0.132 
Imd 40.69 3.822 2945 49.86 | 0.646 0.627 | 0.000 0.000 0.019 
or 25.00 2.268 20.95 31.18 | 0.612 0.595 | 0.012 0.000 0.796 
hop 42.87 3.379 34.91 49.87 | 0.725 0.700 | 0.155 0.000 0.501 
ir 50.34 5.895 37.67 63.73 | 0.804 0.800 | 0.000 0.000 0.304 
ist 40.53 3.142 3280 46.19 | 0.648 0.626 | 0.130 0.000 0.090 
hco 59.50 3.694 50.57 68.30 | 0.703 0.689 | 0.009 0.000 0.845 
g Е 
: 5s. 6 @ § Fb a 8 
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Fig. 2. Dendro- 
30 gramme result- 
ing from cluster 
analysis of mor- 
phometric char- 
acters in round 
40 goby males and 
females of dif- 
ferent SL from 
north-western 
Sea of Azov us- 
50 ing complete 
linkage method. 
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Table 3. Morphometric characters of round goby males (including largest specimens over 15 cm 
long) from north-western Sea of Azov and results of correlation and variance analyses (ANOVA and 
Kruskal-Wallis) to assess dependence on standard length factor. 











Males (n - 65) Correlation analysis о 
Character 
Mean SD Min Max T R ANOVA d 

SL, cm 12.10 3.476 5.4 18.0 — - - - 
?6 SL 

H 23.10 1.436 19.38 27.97 0.408 0.326 0.000 0.004 
h 12.01 0.751 10.45 13.85 0.588 0.602 0.000 0.000 
iH 18.71 1.247 15.85 21.21 0.513 0.443 0.000 0.001 
ih 355.20 0.464 4.40 6.33 -0.251 -0.157 0.008 0.024 
aD 32.61 1.065 29.38 35.61 -0.015 -0.084 0.017 0.072 
pD 19.56 1.181 16.76 22.22 -0.381 -0.370 0.019 0.047 
aP 29.74 1.072 26.82 31.73 0.458 0.401 0.005 0.035 
aV 30.44 1.267 27.30 32.96 0.421 0.405 0.012 0.018 
aA 55.93 1.705 51.93 59.74 0.203 0.189 0.766 0.776 
pow-y 1116.84 1.379 13.39 19.41 0.689 0.617 0.000 0.000 
Роу 5.90 0.570 4.39 7.11 0.434 0.435 0.002 0.006 
V-A 27.13 1.849 22.04 32.31 0.049 -0.018 0.727 0.806 
р1 20.87 1.398 17.86 24.80 -0.292 -0.260 0.131 0.153 
ИО 17.15 1.303 13.54 20.07 0.419 0.386 0.000 0.002 
мр 15.84 0.977 13.99 17.72 -0.373 -0.387 0.005 0.007 
12р 33.02 1.583 29.94 37.69 0.328 0.357 0.248 0.117 
h2D 15.55 0.900 13.29 18.05 -0.239 | -0.253 0.015 0.008 
lA 25.69 1.354 23.19 28.97 0.420 0.426 0.052 0.071 
hA 311.20 0.859 9.71 14.07 -0.631 -0.628 0.000 0.000 
IP 30.18 1.618 26.23 34.33 -0.268 -0.267 0.075 0.076 
iP 312.60 0.920 10.45 14.41 0.810 0.812 0.000 0.000 
IV 23.30 1.285 20.83 27.22 -0.689 -0.659 0.000 0.000 
iV 7.26 0.481 6.06 8.20 0.665 0.636 0.000 0.000 
IC 21.39 1.128 18.89 24.03 -0.348 -0.319 0.094 0.129 
НІ. 27.86 1.021 25.66 30.19 0.424 0.455 0.002 0.003 
% НІ 

һс2 82.49 4.364 69.94 90.83 0.745 0.723 0.000 0.000 
ic 3194.30 6.465 76.69 104.26 0.797 0.779 0.000 0.000 
ao 431.84 2.270 23.93 34.92 0.659 0.662 0.000 0.000 
о 19.97 2.463 14.26 24.15 -0.881 -0.886 0.000 0.000 
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Table 3. Continued. 























Males (n = 65) Correlation analysis E 
Character 
Mean SD Min Max r R ANOVA т 
ор #51.77 1.651 46.01 54.79 0.292 0.305 0.134 0.273 
10 18.88 2.672 10.65 24.55 0.862 0.860 0.000 0.000 
Im 333.84 3.376 24.54 38.25 0.840 0.843 0.000 0.000 
Ња 43.68 4.115 34.46 51.19 0.807 0.816 0.000 0.000 
ог 26.43 2.409 21.26 31.18 0.675 0.667 0.000 0.000 
Бор 45.39 4.171 35.80 53.36 0.866 0.867 0.000 0.000 
ir 1156.25 6.670 41.72 67.07 0.939 0.931 0.000 0.000 
ist 41.12 2.729 34.32 47.09 0.512 0.484 0.000 0.003 
hco 61.64 3.696 50.92 68.87 0.704 0.689 0.000 0.000 
DISCUSSION ing behaviour—while constructing the nest, 


As shown above, the characters that are 
strongly correlated (2 0.75) with SL are not 
the same in males and females. Thus, the 
females have five strongly size-dependent 
characters describing the general body 
shape and the shape of the fins (maximum 
body depth, maximum body width, pecto- 
ral-pelvic distances 1 and 2, and caudal fin 
length) while the males only one (width of 
pectoral fin base). It can be hypothesised 
that such a correlation between the width 
of the pectoral-fin base and SL in males is 
caused by their spawning behaviour — a 
male, guarding the nest, actively aerates the 
eggs moving the pectoral fins. 

The head parameters strongly correlated 
with SL are more numerous in the males 
(seven vs. five in the females). Three char- 
acters are shared by the males and females 
(maximum head width, interorbital dis- 
tance, and mouth width); two characters 
are specific for the males (horizontal eye 
diameter, cheek depth, upper jaw length, 
and lower jaw length) and two characters 
are found only in the females (head depth 
at nape and isthmus width). It is probable 
that a strong correlation between the cheek 
depth and length of jaws in males and their 
body length, are also relevant to the spawn- 


the male carefully cleans it from the "gar- 
bage" moving it outside the nest with the aid 
of the mouth. Later, when guarding the laid 
eggs, the male scares other fish demonstrat- 
ing aggressiveness by inflating the cheeks. 

The statistical methods used in the 
study proved to be appropriate for discov- 
ering variability of different nature. 

The applicability of methods of paramet- 
ric statistics is determined primarily by the 
normal distribution of the random variable 
distribution though most of actual distribu- 
tions of measurements of biological objects 
are not normal (Orlov, 2002). So, results of 
studies using parametric statistics, especial- 
ly the frequently used t-test, for non-normal 
distributions and small compared samples, 
should be accepted critically. 

From the theoretical point, thousands 
of observations are needed for any reliable 
verification of the normality. However, if 
the homogeneity of two samples is to be 
studied using a t-test (with an a priori as- 
sumption of the equality of variances), then 
the effect of departures from the normality 
decreases with the increase of the sample 
size (Orlov, 2002). In practice, the concept 
of “large” and “small” samples are used in bi- 
ological and fisheries research (Aksyutina, 
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1968; Lakin, 1980), a large enough sample 
containing over 30 observations when, with 
certain reservations, parametric statistics 
can be applied. 

The methods of parametric and nonpar- 
ametric statistics used in the present study 
resulted in very similar estimates that is ap- 
parently connected with a sufficiently large 
size of the analysed samples. In general, 
when the normal character distribution is 
considered, the parametric tests are more 
powerful than the nonparametric ones. So, 
in all cases when compared samples are tak- 
en from normally distributed populations, it 
could be recommended to prefer the para- 
metric tests. In case of considerable depar- 
tures from the normality, nonparametric 
tests should be used. 

For the analysis of size-dependent vari- 
ability, as shown by our results, it was not 
sufficient to base conclusions on correla- 
tion test statistics only. A comparison of 
the data based on the analysis of variance 
with the results of the correlation analysis 
revealed certain differences. Even in cases 
of a statistically significant correlation, 
but at relatively not very high values of the 
correlation coefficient (less than 0.5), the 
analysis of variance, in many cases, does not 
confirm the existence of such a relationship 
(Tables 2, 3). In other words, a statistically 
significant but weak correlation revealed in 
the data set for each individual character, is 
not necessarily evident from the in-group 
analysis of the character. For example, the 
results of the correlation analysis of size- 
dependent variability in round goby males 
5-18 cm long showed the statistically sig- 
nificant correlation for caudal peduncle 
length (r = —0.292), length of second dorsal 
fin base (r — 0.328), length of anal fin base 
(т = 0.420), pectoral fin length (r = —0.268), 
caudal fin length (r ^ —0.348), and postor- 
bital distance (r — 0.292) but the results 
of the ANOVA demonstrate the absence of 
correlation of these morphometric charac- 
ters with SL. 

This clearly supports the known idea 
(e.g., Aivazyan et al., 1985) that the corre- 
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lation analysis though allows to reveal the 
presence and strength of statistical rela- 
tionships between two variables but is not 
sufficient to prove the existence of a causal 
link between characters (SL and relative 
morphometric characters in our case). A 
more reliable conclusion can be only ob- 
tained when a variance analysis is also ap- 
plied because it estimates in-group variabil- 
ity of the factors and their joint effect. This 
is especially true for ambiguous results, for 
example, 1) if results of the used statistical 
tests (e.g. U-test and t-test) do not coincide 
as in the case of maximum body depth and 
pectoral-pelvic distance 2 when the females 
and males were compared (Table 2), or 
2) analysing correlation between variables 
(characters), especially in the range of low 
values of the correlation coefficient. 

It should be noted that in recent years, 
some authors (Kováé, Syriová, 2005; 
Eavrincíková et al., 2005) successfully ap- 
plied a regression analysis using geometric 
morphometric methods for the study of de- 
velopmental variability in gobies, including 
the round goby. 

Summarising conclusions on the meth- 
odology used for a morphometric analysis 
are as follows. 1. In order to reach a cor- 
rect assessment, a study of morphological 
variability, especially of sex- and size-de- 
pendent differences, should use appropri- 
ate methods selected with a particular at- 
tention to the actual characteristics of data 
to be analysed (distribution of variables, 
sample size). 2. Evaluation of correlation 
should be based not only on values of a 
correlation coefficient calculated for the 
entire set of data per character but, also, 
on results of an in-group analyses of the 
pattern and extent of correlation within 
the range of data proper, e.g., within size 
groups or classes. This approach is most 
fully realised in the methods of variance 
analyses. The approach of searching for the 
most informative size range could be useful 
for further studies of infraspecific variation 
and comparison between morphologically 
close gobiin species. 
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